PSEUDOPROCESSES ON A CIRCLE AND RELATED POISSON 

KERNELS 



ENZO ORSINGHER AND BRUNO TOALDO 



Abstract. Pseudoprocesses, constructed by means of the solutions of higher- 
order heat- type equations have been developed by several authors and many re- 
lated functionals have been analyzed by means of the Feynman-Kac functional 
or by means of the Spitzer identity. We here examine pseudoprocesses wrapped 
up on circles and derive their explicit signed density measures. We develop the 
Fourier series of their laws in the case of even-order and odd-order pseudopro- 
cesses separately. We observe that circular even-order pseudoprocesses differ 
substantially from pseudoprocesses on the line because - for t > t > 0, where 
f is a suitable n-dependent time value - they become real random processes. 
By composing the circular pseudoprocesses with positively-skewed stable pro- 
cesses we arrive at genuine circular processes whose distribution, in the form 
of Poisson kernels, is obtained. The distribution of circular even-order pseudo- 
processes is similar to the Von Mises (or Fisher) circular normal and therefore 
to the wrapped up law of Brownian motion. The last section of the paper is 
related to the circular Fresnel pseudoprocess and we examine its composition 
with stable subordinators. This leads to combinations of Poisson kernels. 



1. Introduction and preliminaries 

Pseudoprocesses are connected with the fundamental solution of heat-type equa- 
tions of the form 

d 9" 

— u„(a;,t) = c„— — u„(x,t), a; G M, t > 0, n G N, (1.1) 
ot ox"- 

where 

{(— l)t+i, for even values of n 

(1.2) 
±1, for odd values of n, 

subject to the initial condition 

m(x,0) S{x). (1.3) 

For n > 2 the fundamental solutions to (1.1) are sign-varying. By means of a 
Wiener-type approach some authors (see for example Albeverio et al. [1], Daletsky 
[10], Daletsky and Fomin [11], Krylov [16], Ladohin [19]) have constructed pseudo- 
processes which we denote by X{t), t > oi X„(i), if we specify the order of the 
governing equation. In these papers the set of real functions x : t G [0, cxi) — >■ x(t) 
(sample paths) and the cylinders 

C = {x{t) : aj < x{tj) <bj, j = 1, - ■ ■ ,n} (1.4) 
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have been considered. By using the solutions u„ to (1.1) the measure of cylinders 
is given as 

/•bi /•&„ " 

^J■n{C) = / dxi--- dx„ Yl '^n{Xj ~ Xj^i,t^ - tj^i) . (1.5) 

ai •'On j — i 

In (1.5) we denote by u„ 

un{x,t) = — / d^e-'«^e^"(-*«)"*, (1.6) 

in case the integral above converges. The measure (1.5) is extended to the field 
generated by cylinders (1.4) for fixed ti < ■ ■ ■ < tj < ■ ■ ■ < tn- The signed measure 
obtained in this way is Markovian in the sense that 

{X {t + T)e 'b\^t} = Mx(T) {Xit) e ®} , (1.7) 

where ^Ft is the field generated as 

Tt = (T{X{h) e ,X(t„) e ®„}. (1.8) 

More information on properties of pseudoprocesses can be found in Cammarota 
and Lachal [9], Lachal [17] and Nishioka [22]. For pseudoprocesses with drift the 
reader can consult Lachal [18]. 

In this paper we consider pseudoprocesses on the ring ^ of radius one, denoted 
by 8(t), t > 0, whose signed density measures are governed by 

(§-^Vn{e,t) ^ c„^Vn{e,t), [0,27r),t>0,n>2, 

\vn{e,o) = s{e). 

The signed measures of pseudoprocesses on the line X{t), t > 0, and those on the 
unit-radius ring, 0(i), t > 0, can be related by 

oo 

{Q{t)ede}^ [j {X{t) £d{e + 27mT)}, 0<6'<27r. (1.10) 

This means that the pseudoprocess O has sample paths which are obtained from 
those of X by wrapping them up around the circumference Counterclockwise 
moving sample paths of Q correspond to increasing sample paths of X. 

For n = 2 we have in particular the circular Brownian motion studied by Hart- 
man and Watson [14], Roberts and Ursell [27], Stephens [28]. The pseudoprocesses 
running on are called circular pseudoprocesses and are denoted either by 6(i), 
t > 0, or On{t) if we want to clarify the order of the equation governing their distri- 
bution. We analyze separately the even-order case and the odd-order one because 
they pose qualitatively different problems. In view of (1.10) we can write 

oo 

V2n{0,t) = J2 U2n{0 + 2mTT,t), < 9 < 2tT . (1.11) 

m— — oo 

Equation (1.11) shows that the solution to (1.9) can be obtained by wrapping up 
the solution to (1.1) which reads 

1 2 

U2n{x,t) = — d^e-^^'^e-^'"*, xeM,<>0. (1.12) 
27r J_oo 

The function (1.12) has been investigated in special cases by Hochberg [15], Krylov 
[16], Nishioka [22] and more in general by Lachal [17; 18]. The sign- varying structure 
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of (1-12) has been discovered in special cases by Bernstein [7], Levy [20], Polya [26], 
as early as at the beginning of the Twentieth century and has been more recently 
studied also by Li and Wong [21]. 

The Fourier series of (1-11) has the remarkably simple form 

1 1 2 

V2n{0,t) = — + e-''^"*coske, 6le[0,27r). (L13) 

k = l 

For n = 1 we obtain the Fourier series of the law of the circular Brownian motion 
(see [14]). The function 

1 1 °° 

V2{9,t) ^ — + e"'"'' cos k9 (1.14) 

k=l 

is similar to the Von Mises circular normal 
where 

/ T* \ 27+m 1 

W-)-g(l) .!r(„. + ,-i) 

is the m-th order Bessel function. The relationship between (1.14) and (1-15) is 
investigated in the paper by Hartman and Watson [14]. The Von Mises circular 
normal represents the hitting distribution of the circumference !^ of a Brownian 
motion with drift starting from the center of The planar Brownian motion 
{R{t), 4'(i)), t > 0, with drift k = (fci,A:2), ||fc|| = k, has transition function 



Pr{i?(i) e dp, *(<) e dip} = -^e-'^e-'^eP'"'°'''^dpdip (1.17) 

27ri 

and marginal 



Therefore 



Pv{R{t) e dp} = ^e-'^e-^Io{pk)dp. (1.18) 

^pk cos if 

Pr {*(t) € dip\Rit) e dp} = ^^^^^^ dip (1.19) 

and for p = 1 coincides with (1.15). 

The analysis of the pictures of V2n{d, t) for different values of t and different val- 
ues of the order 2n, n S N, shows that the distributions (1.13) after a certain time 
become non-negative. This means that pseudoprocesses on the circle "R^ behave dif- 
ferently from their counterparts on the line and rapidly become genuine processes. 
Furthermore we remark that in small initial intervals of time the circular pseudo- 
processes have signed-valued distributions with a number of minima which rapidly 
unify into a single minumum (located at = vr) which for increasing t upcrosses the 
zero level. This is due to the fact that in a small initial interval of time the effect 
produced by the central bell of the distribution has not yet spread on the whole 
ring 01. 

The value of the absolute minimum of V2n{0, t) for t > t has the form 

1 1 °° 

^^2„(^,i) = + (1-20) 

ZTT TT ^ — ' 

k=l 
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The graph of functions V2n{d,t) shghtly differ from that of the density of circular 
Brownian motion as shown in figures 1 and 2. The term fc = 1 in (1.13) is the 
leading term of the series and the form of the distribution f2n(^, t) is very close to 
that of + e~* cos 9. 



Figure 1. The distributions of the fourth-order circular pseudo- 
process for different values of t 





Figure 2. The distributions (for t = 1) of the circular pseudopro- 
cesses of various order 2n 
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The odd-order case is much more complicated because the solutions to equation 
(1.1) are asymmetric (with asymmetry decreasing for increasing values of the order 
n). Some properties of solutions to odd-order heat-type equations can be found in 
[17; 18]. In the present paper the wrapped up solution to (1.1) gives the fundamental 
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solution of (1.9) as 

oo 

V2n+l{0,t) = ^ M2„+l(0 + 2fc^,i), 0€[O,27t), 

whose Fourier series reads 

V2n+liO,t) 



k — — oo 



1 1 



2tt 



-^cos(fc2"+4 + fc6'). 



(1.21) 



(1.22) 



fe=i 



We note that for n = 1 the series (1.22) becomes a discrete version of the solution 
to (1.1) which reads 



U3{x,t) 



where 



Ai(x) 



1 










" Wt 






/•oo 












/ COS 






Jo 







da 



(1.23) 
(1.24) 



is the Airy function. 

We consider also the wrapped up stable processes 6^'^(t), t > 0, and the re- 
lated governing space-fractional equation. In particular we show that the law of 
Q"^!^ (2~^t), t > 0, is the fundamental solution of the space-fractional equations 



e [0,27r) , i > 0,/3 e (0, 1] 



and has Fourier expansion 



1 

2tt 



1 



oo / 2\t- 



cos mb 



(1.25) 



(1.26) 



The fractional operator appearing in (1.25) is the one-dimensional fractional Lapla- 
cian which can be defined by means of the Bochner representation (see, for example, 
Balakrishnan [3], Bochner [8]) 

We show that formula (1.26) coincides with the distribution of the subordinated 
Brownian motion on the circle, *B (i?^(t)), t > 0, where H^{t), i > 0, is a stable 
subordinator of order f3 G (0,1] (see, for example, Baeumcr and Mecrschacrt [2]). 
Furthermore we notice that 

<B (2 i?^ (<))'= 62^(t), t>0, (1.28) 

©^^(t), t > 0, is a symmetric process on the ring with distribution which can be 
obtained by its symmetric stable counterpart on the line as 



oo 



^i^(e+27n7r)^-t\i\^f 



1 

2^ 



cos k9 



k=l 



(1.29) 
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For /3 = i we extract from (1.26) the Poisson kernel 

vl{e,t) = , . (1.30) 

ivr 1 4_ g-i\/2 _ 2e v^cos6' 

The composition of the circular pseudoprocesses Qn(t), t > 0, with positively- 
skewed stable processes of order ^, say H^{t), t > 0, leads also to the Poisson 
kernel. In particular, we show that 

Q2n(H^-{t))^de\ = - ^ 0e[O,27r). (1.31) 

V / J 27r 1 + e — 2e 'cosfc' 

In the odd-order case the result is different, depends on n and has the following 

form for 6* € [0, 27r) 

{e2„+i (if5T^(t)) ed^^} = 

, (1.32) 

2^ 1 + e"^*'=°' 5(5^^+1) - 2e"*™"5(2;^+T) cos (6> + tsin 2(2,^^1) ) 

The composition of pseudoprocesses with stable processes therefore produces gen- 
uine r.v.'s on the ring as it happens on the line (see Orsingher and D'Ovidio 
[24]). We note that the distribution of the composition in the even order case is 
independent from n (formula (1.31)), while in the odd-order case the Poisson kernel 
obtained depends on n and has a rather complicated structure. For n ^ 00 the 
kernel (1.32) converges pointwise to (1-31) since the asymmetry of the fundamental 
solutions of (1.1) (as well as that of their wrapped up counterparts) decreases. The 
result (1.31) offers an interesting interpretation. The Poisson kernel (1.31) can be 
viewed as the probability that a planar Brownian motion starting from the point 
with polar coordinates (e~*, 0) hits the circumference ^ in the point (1, Q) (see Fig. 
4a) . Therefore this distribution coincides with the law of an even-order pseudopro- 
cess running on the circumference and stopped at time H^{t), t > 0. This result 
is independent from n and therefore is valid also for Brownian motion. A similar 
interpretation holds also for circular odd-order pseudoprocesses taken at the time 
fj 2,1+1 (i)^ t > 0, but starting from the point with polar coordinates (e"""*, 6„i), 
where a„ = cos7r/(2(2ri -t- 1)) and 6„ = sin7r/(2(2n + 1)). 

A further generalization of the above framework is obtained by considering the 
fractional higher-order diffusion equations (studied on the whole real line by Beghin 
[4]) 



§^v-2n{0,t) = C2n^v-^^{e,t), G [0, 27r) , t > 0, n € N, 



(1.33) 



The solution to (1.33) can be written as 



^ ^1 + 2 £ E^^i (-rm^") cos me^ , (1.34) 
where 

L^it) ^ M{s>0: H''{s)>t}, t > 0, (1.35) 

is the inverse of the stable subordinator H^{t), t > 0. The fractional derivative 
in (1.33) must be meant in the Dzcrbayshan-Caputo sense and E^ i represents the 
Mittag-Leffler function. Clearly for i/ = 1 we retrive result (1.13). 
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The Fresnel pseudoprocess F{t), t > 0, related to the equation of vibrations 
of rods was introduced and investigated in Orsingher and D'Ovidio [23]. We here 
study the corresponding version on the ring say ^(t), t > 0, and obtain the 
fourth-order Poisson kernel representing the distribution of 5^ (rt), t > 0, where 
Tf = inf {s > ; B{s) = t}. Circular Fresnel pseudoprocesses related to fractional 
equations of the form 

n^.t) ^ --r—r{0,t), 0e[O,27r),t>O, z.e(0,l], (1.36) 



df2^' y ' / 22 9614 
are analyzed and the corresponding densities 



1 1 / i2i/y^4\ 



2tT TT^ \ 22 

k=l ^ 

are obtained. The composition of the circular Fresnel pseudoprocess with Brownian 
first passage times leads to the superposition of Poisson kernels. 

2. Pseudoprocesses on a ring 

In this section we consider pseudoprocesses Q{t), t > 0, on the unit-radius cir- 
cumference ^, whose density function Vn{0,t), 9 £ [0,2n), t > 0, is governed by 
the higher order heat-type equation 

i§-^vn{0,t) ^ Cn-§^v^{e,t), 9e[0,2TT),t>0,n>2, 
\vn{9,0) = 6{9). 

The pseudoprocesses 0„ have sample paths obtained by wrapping up the trajec- 
tories of pseudoprocesses on the line M. Increasing sample paths on !^ correspond 
to counterclockwise moving motions on the ring The structure of sample paths 
of pseudoprocesses has not been investigated in detail although some results by 
Lachal (Theorem 5.2, [18]) show that there is a sort of "slight" discontinuity in 
their behaviour (this is confirmed by Hochberg [15]) and the fact that the reflection 
principle fails (Beghin et al. [5], Lachal [17]). 

It must be considered that the wrapping up of paths and the corresponding 
density measures produces in the long run genuine processes (with non-negative 
measure densities in the case n is even). Our first result concerns the distribution 
of e„(i), t > 0. 

Theorem 2.1. The solutions to the n-th order heat-type equations (2.1) reads 

iv2ni9,t) = 2^ + i^^^^e-''^'"*cosfc0, for C2n = n > 1, 

\v2n+i{9,t) = 2^ + i^^^^cos(fc2»+ii + fc0), /orc2„+i = (-1)", n> 1. 

(2.2) 

Proof. We can obtain the result (2.2) in two different ways. We start by considering 
the even-order case where the wrapping up of the solutions to (1.1) leads to 

V2n{9,t) = U2n{9 + 2mTT,t) = ^ / e-'('^+2'"-)«e-«"*. (2.3) 

m— — oo m— — oo 

The Fourier series expansion of the symmetric fmiction V2n{p^t) has coefficients 



1 

ttk ^ — d9 cos k9 

Jo 



U2n{9 + 2rmT,t) 
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= 2 E 



oo /-m+l 



1 
TT 
1 

2^ 
1 

2^ 



— OO 
CX3 



cos/cz 



— CX3 
OO 



1 

2^ 



— OO 
OO 



d^e-^ * [(5(^ - fc) + + A:)] = 



(2.4) 



This shows that the first result of (2.2) holds. An alternative derivation of V2,i{d,t) 
is based on the method of separation of variables. Thus under the assumption that 
V2ni0,t)=T{t)TP{e) we get 



r(i) v(^?) ^ ' ^ ' 



(2.5) 



In order to have periodic solutions we must take integer values of /3 and thus the 
general solution to (2.1) becomes 



It) ^ ^ AkC''^"* COS k9 = + e^''^"^ COS ke. (2.6) 



fc — — OO 

The initial condition 



fc=i 



1 °° 

f2„(e,0) = 5 {6) = +- Vcosfce 



(2.7) 



implies that Ak = ^ . In the odd-order case, by proceeding as in the previous case 



we have that the Fourier coefficients of 

OO 

m— — OO 

become 



(2., 



ak 



1 

2^ 
1 

2tt 
1 

2^ 



1 

2^ 



['5(e -k) + + k)] 



TT 



(2.9) 



In a similar way we have that 

.(-i)"(-«?)'"+'t 



bk = 



27T 
1 

27ri 



OO 
OO 



1 



dz 



27ri 

dCe-'^""'"'['5(C-fc)-'5(e + fc) 
and thus the expression (2.2) oi V2n+i{0,t) follows. 



■sinfc2«+ii, 



(2.10) 

□ 
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In a recent paper by Orsinghcr and D'Ovidio [24] the authors have shown that 
the solutions to (1.1) have the following probabilistic representation 

(u2n{x,t) = ±^E{sm{xG^" (j))}, forc2„ = (-l)"+i,n>l, 

\u2n+i{x,t) = ^]E{e-''"^^"^'(7)sin(a„:EG2"+i(l))}, for C2„+i - (-1)", n > 1, 



where C (t ^) is a generalized gamma r.v. with density 



X > 0, < > 0, 7 > 0, 



and 



(2.11) 
(2.12) 

" 2(2n + 1) ' " 2(2n + 1) " ^^'^^^ 

We now confirm the results of Theorem 2.1 by means of the wrapped up versions 
of the probabilistic representations (2.11) which become 

• t) = E™=-oo ^F(«4k^E{sin((0 + 2™.)G2"(i))}, 



TT 



V2v 



V2n+l{e,t) = J2'^= 



m=— oo ir(0+2m7r) 



E 



|g-6„(e+2™.)G-+^(i) ^^^(g _^ 2m^)G2"+i {}))} 

(2.14) 



We prove that the Fourier series expansion of (2.14) coincides with (2.2). We 
start first with the even-order pseudoprocesses for which only the cosine Fourier 
coefhcients Ofc are necessary. 



Proof of (2.15) in the symmetric case. We have that 



OO 



m— — CX3 



de 



^2 



9 + 2mTT 



COsfcz . / 2ri 

dz sm zG 



E <^ sm{0 + 2mTT)G 



2n 



COS k9 



— oo 
oo 



dz 



z 

sin 



in(fcz + zG2" (i)) sin(zG2" (i) - fcz) 



-e(i 

TT L 



-G-(i)<fe<G-(i)] 

where we used the fact that 

sin ax 

dx 



> k 



- J 2' 



if a > 0, 
if a < 0. 



(2.15) 



(2.16) 



The calculation (2.15) shows that the density of even-order circular pseudoprocesses 
can be viewed as the superposition of sinusoidal waves whose amplitude corresponds 
to the tails of a Weibull distribution. □ 

Proof of (2.14) in the asymmetric case. In the odd-order case we need both the 
sine and cosine coefficients of the Fourier expansion because the signed laws are 
asymmetric. For 



V2n+l{d,t) = ^ 



1 



m— — OO 



2m7r) 



2m7r)G2"+i 
(2.17) 
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the Fourier coefRcients become 



\ak = icOsfc2n+lt, 

jfofc = -isin/fc2"+it 
We give with some details the evaluation of (2.18) 

Sin (a„zG-+i (l))} 



(2.18) 



^ — ^ Q -\- 2m7T 

m— — oo 

1 ^ / ^ cosfcz 
E / dz — e 

1 

2^^2 
1 



2rmr)G^'' 



2n+l 



oo 
oo 



sm 



2z27r2 



E 



____ (z (a„G^"+i (i) + fc)) + sin (z (a„G^"+i (f) - fc)) 

z 

g..(a„G- + i(i)+fc)-6„.G- + i(i) „ ^-,.(a„G-+i(i)+fc)-b„.G-+i(i) 



Z 



oo 



^g«(a„G-+^(i)-fc)-6„.G- + i(i) _ g-..(a„G-+^(i)-fc)-6„.G- + i(i) 

(2.19) 

By considering the following integral representation of the Heaviside function 



ny{x) = - 



1 

2^ 



(2.20) 



the coefficients in (2.19) become 



{2n + l)t 



dw u;2"e-*'^'"^' [•Hfe(w(a„ - ib„)) - ■Hfc(-w(a„ + ib,,)) 



2^ JO 

+Hk{w{an + ibn)) - 'Hfe(-w(a„ - i6„)) 



i(2n+ l)i 
2^? 



H,(u;) + ^i^I^±ll^ / dww'-e^-^"^''nk{^w) 



2tt 



i{2n + l)t 



2tt 

i{2n + l)t 
2^ 



dw \ 



dww''''e-''"'"^''Hk{w) 



i{2n + l)t 
2tt 



27r v / 



-cos fc2"+lt. 



(2.21) 



In order to justify the last step we can either take the Laplace transform with 
respect to t (see for example Orsingher and D'Ovidio [24]) or we can apply the 
following trick 

C^O 27T VA Jk 

(2.22) 

The coefficients bk (2.18) can be obtained by performing similar calculation. □ 
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2.1. Circular Brownian motion. The circular Brownian motion < > 0, 

has been analyzed by Roberts and Ursell [27], Stephens [28] and also by Hartman 
and Watson [14]. In a certain sense it can be viewed as a special case of symmet- 
ric pseudoprocesses on the ring The distribution of *B(t), t > 0, has Fourier 
representation 




l + 2^e-Vcosfc6l , 6le[0,27r), (2.23) 
fc=i / 

and can be also regarded as the wrapped up distribution of the standard Brownian 
motion 

1 °° 2 

p^{e,t) - ^ E e-""^^. (2.24) 

^ m— — oo 

Formula (2.23) corresponds to n = 1 of (2.2) for the even-order case with a suitable 
adjustement of the time scale. The law (2.23) can be obtained directly by solving 
the Cauchy problem 

^^tP's{o,t) = i-^p^{0,t), 0e [o,2^),i>o, 

\p<s{0,O) = 6(0). ^ ■ ^ 

or as the limit of a circular random walk as in [28]. The distribution of the circular 
Brownian motion is depicted in Figure 2 and looks like the Von Mises circular 
normal (this is the inspiring idea of the paper by Hartman and Watson [14] in 
which the connection between the two distributions is investigated). For t ^ oo 
the distribution of *8(t), i > 0, tends to the uniform law. 
We note that 

oo _ (2fc+l)^t 

fe=0 

and therefore 

Pr(-^ < «8(t) < ^1 < ^ + ^e"^ valid for t > -21og^ w 0.209. (2.27) 
12 2 J 2 TT 4 

The relationship between circular Brownian motion *B(t), t > 0, and Brownian 
motion on the line B{t), t > 0, 

oo 

{^{t)ede}= [j {B{t)ed{e + 2rmr)}, 6* € [0, 27r) , (2.28) 

m— — oo 

permits us to derive the distribution of 

max |«B(i)| , t > 0, (2.29) 

0<s<t 

that is the distribution of the maximal distance reached by the circular Brownian 
motion from the starting point. Of course the sample paths overcoming the angular 
distance n at least once are assigned tt as maximal distance which therefore has a 
positive probability (converging to 1 as time tends to infinity). 

Proposition 2.2. For the maximal distance (2.29) we have that 
Pr I max |Q3(s)| < e\ = [ Pr I-^ < min B(s) < max B(s) < 

|^0<S<t J J_g I 0<S<t 0<S<t 
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{y-4mg)^ OO ( - « + 20 ( 2m - 1) ) ■= 

e 2t e 2t 



\m— — oo ^ m— — CX3 

. (l-'2r)0 



(-l)"^ / dw^^. (2.30) 

^ ^ / (l+2r)0 ^/2^ ^ ^ 



The related first passage time of circular Brownian motion has density which has 
the following form 

Pr {Te edt} - ^ Pr | max |»(s)| < 6* 1 dt 
dt o<s<t 



^ ^ , 9{l - 2r) + ^ ^ , e{l + 2r) 



Curiously enough the factor 6e~^ /y2'Kt'^ coincides with the first passage time 
through of a Brownian motion on the hne. 

3. Fractional equations on the ring 3?, and the related processes 

In this section we consider various types of processes on the unit radius circum- 
ference 

3.1. Higher-order time- fractional equations. We start by analyzing the pro- 
cesses related to the solutions of time-fractional higher-order heat-type equations. 
We consider the time-changed pseudoprocesses 92„ {L'^{t)), t > 0, where 

L^it) = inf {s > : H^is) > t} (3.1) 

and where H''{t), t > 0, is a positively skewed stable process of order i' £ (0, 1]. 
We notice that the Laplace transform of the distribution l,y{x,t) of (3.1) reads (see 
for example Orsingher and Toaldo [25]) 



dxe-^''l^{x,t) = E^^i{--ft'') (3.2) 



where 



^ r(i^j -I- 1) 

is the Mittag-Leffler function. For pseudoprocesses related to time-fractional equa- 
tions we have the next theorem. 



Theorem 3.1. The solution to the problem, for v e (0, 1], n G N, 

fi^«2J^,i) = -(-^i^)"«2n(^,i), [0,27r),i>0, 

U2^„(0,o) - m. 

is the univariate (signed) distribution o/ 02„ {L'^{t)), t > 0, which reads 



(3.4) 



if °° 
\ fc=i 



V2niO,t) = — 1 + 2 yE,,, (-k'^f') coske . (3.5) 
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The time-fractional derivative in (3.4) must be understood in the Caputo sense, that 
is 



1 



zny , , Y(\-v) Jo {t-sy 
Proof. The law of Oa™ (L^it)), i > 0, is given by 



0<iy <l. 



(3.6) 



t) = 



^ + 2 £ ^^,1 {-k^^t") cos k0^ 



cos k9 ly{s,t) 



(3.7) 



Since, Vfc > 1, we have that 



9" 



dt^ 



n2n 

(-fc2"r) cos/fc^ = -k^^'E.^i {-k^'H"") coske = (-l)"+'^^^.,l (-fc'"t") cosfc^, 



and therefore we conclude that (3.7) satisfies the fractional equation (3.4). 



(3.8) 

□ 



Remark 3.2. For n ^ I, formula (3.7) becomes the distribution of subordinated 
Brownian motion S {L'^{t)), t > Q. For v = 1 we retrieve from (3.5) the solutions 
(2.2) of the even-order heat-type equations on 

3.2. Space-fractional equations and wrapped up stable processes. The fol- 
lowing Theorem represents the counterpart on !^ of the Riesz statement on the 
relationship between space-fractional equations and symmetric stable laws (for the 
non-symmetric case see the paper by Feller [12]). 

Theorem 3.3. The law of the process 03 (i?'^(<)), t>Q, is given by 



Pi{0,t) 



1 

2^ 



1+2^:, 



' cos kO 



(3.9) 



and solves the space-fractional equation, for /3 G (0, 1], 



jkp'^i9,t) = -[-^^) pi{e,t), 

,{0,O) = S{9). 

The fractional one- dimensional Laplacian in (3.10) is defined in (1.27) and H^{t), 
t > 0, is a stable subordinator of order (5 G (0, 1] . 

Proof The law of <B {H'^{t)), t > 0, is given by 



e [0,27r), t > 



(3.10) 



piio,t) 



dsp<si0,s) hp{s,t) 



1 

271^ 



1 



;2_]e V 2 J 



' COS kO 



k=l 



(3.11) 



where p<s is the law of circular Brownian motion and hp is the density of a positively 
skewed random process of order /3 S (0, 1]. In order to check that (3.9) solves (3.10) 
it is convenient to write the fractional derivative appearing in (3.10) as 

P ■ n ^oo / / -, o9 X X -1 



19^ 
2d(P 



sin 7r/3 



A + 



2d(P 



\P dX 



r(/?)r(i-/3) 7o 



e-"^-<---^) dudX 
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w-'^-ie~"v du. (3.12) 



r(-/3) Jo 

From (3.12) we have therefore that 

cosfc6> = / duu-f^-^e cosA:6i 

2 oO J r (— p) Jo 

-r(-/3)X j^, ji [ 2) 902.^°^^' 

= — — — / duu ^ y , cos kO 



1 



r(-/3) 



2 , cosfc6i, (3.13) 

and this shows that (3.9) satisfies (3.10). □ 
Remark 3.4. Another way to prove that 



coske = — coske (3.14) 



can be traced in the paragraph 4-6, page 428 of Balakrishnan [3], which confirms 
our result. 

Theorem 3.5. For the wrapped up version, say &^{t), t > Q, of the symmetric 
stable processes S'^^{t), t > 0, with characteristic function Ee*^'^ = e^*'^' , we 
have the following equality in distribution 

e^^it) S {2HP{t)) '= S {H^ {2l't)) , t>0. (3.15) 

Proof. The density of 6^^(t), t > 0, must be written as 

Pe..(M) ^ ^^i^^^ ^ ^ £ r.|.~'e-'l^l" ,3.0) 

m — — 00 ^ 

The Fourier expansion of (3.16) becomes 

00 

PQ2i3{9,t) = Y+^afeCosfc6' (3.17) 

k=l 

where 

(.277 1 00 ,.00 



- / dO— J2 / d?e-'«('^+^™'')e-*l«l''cosA:^ 

1 ^ /'oo /•m+l 



m— — oo 



{2^f 



dCe-*l«l / dye 
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= _L / d^e-*l«l^'[<5(e-fc)+^(e + fc)] - -e-*'^^'. (3.18) 

J-oo 

This permits us to conclude that 

Pe2p{9,t) = — + - Ve-'="'*cosA:6'. (3.19) 



27r TT 

A,-l 



□ 



While the integral in (3.16) (representing the Fourier inverse of symmetric stable 
laws) cannot carried out, its circular analogue can be explicitely worked out and 
leads to the Fourier expansion (3.19). 

Corollary 3.6. In view of the results of Theorems 3.1 and 3.3 we have that the 
solution to the space-time fractional equation, for j3 G (0, 1], 

§^p-^^{e,t) = -(-ij^)%^^(0,t), ee[o,2n),t>o 
y/i0,o) ^ s{9). 

can be written as 



(3.20) 



fc=i \ 



Ps'^C^'O - ^ + I - ( V ) ^'Icosfcfl, (3.21) 

and coincides with the law of the process 

f-'-'^{t) ^ ^{H^ {L-'it))) , t>0. (3.22) 

In (3.22) is a stable subordinator of order /3 e (0, 1] and is the inverse of 
as defined in (3.1). 

Proof. Here we only derive the distribution of ^'^'^{t), t > 0. We have that 

/■OO /"OO 

Pr {^"''^t) e d0} =d0 Pr {<B(s) e d9} / Pr {H^{w) G ds} Pr {L^it) e dw} 
Jo Jo 



de 


de 


2^ ' 


^ — 

TT 


de 


de 




f — 

TT 


de 

2tt 


1 + 



fe=l 

OO 



/ Ve"^~^ ""cosfce Pr{L''(i) € dw} 



k=l 

1 + 2 22coskeE^^il - i—j t 



(3.23) 



where in the last step we applied (3.2). □ 

4. From pseudoprocesses to Poisson kernels 

In this section we show that the composition of pseudoprocesses of order n run- 
ning on the circumference !^ with positively skewed stable processes of order ^ leads 
to the Poisson kernel. This is the circular counterpart of the composition of pseu- 
doprocesses with stable subordinators which leads to Cauchy processes. In both 
cases pseudoprocesses stopped at i? " (t), t > 0, yield genuine random variables. 

We distinguish the case where n is even from the case of odd-order pseudopro- 
cesses. We have the first result in Theorem 4.1. 
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Theorem 4.1. The composition 02„ \^H 2n{t)j, t > 0, of the pseudoprocess 
with the stable process H^{t), t > 0, has density 

de 1 - e-2t 
2^ 1 + e-2* - 2e-* cos 9 ' 

and distribution function 



distribution function 

Pr{e2„ (i/^w) <^} = 



n e N, 



(4.1) 



arctan ^f^f^ tan |^ 

\^l-e-* 



arctan 



2^ 

tan ; 



9 € [0, tt] , 
e (7r,27r) 



(4.2) 



which are independent from n. 
Proof. We have that 



Pr{e2„ (H^it)^ e de^ ^'^^ Jo 



1 1 v-^ -^.2 

+ - > e 

A:=l 



27r 



' COS kO 



hi{s,t) 



= d9 

_ dO 
~ 2^ 
d9 



1 

2^ 

14 



cos A:6'e' 



-kt 



k=l 

,-t+ie 



„-2t 



2tt 1 



2e- 



(4.3) 



The result (4.2) is derived by applying formula 2.552(3) page 172 of Gradshteyn 
and Ryzhik [13] 

Va^ — 62 tan | 



a + 6 cos X ^a? - ^2 



arctan 



(4.4) 
□ 



Remark 4.2. T/ie Poisson kernel (4.1) can 6e interpreted as the distribution of the 
process B (Tj^) where B is a planar Brownian motion and%<K^ = inf {t > : B{t) Cz !^}. 
In the case of Theorem 4-1 the planar Brownian motion starts from the point 
(e~*,0). Therefore we have that 



law 



s(TO - e^niH^it) 



t > 0. 



(4.5) 



This means that a pseudoprocess running on the ring !^ and stopped at a stable time 
H^{t), t > Q, has the same distribution of a planar Brownian motion starting 
from (e~*,0) at the first exit time from the unit-radius circle. The result (4.5) 
holds for all n & N and represents the circular counterpart of the composition of 
pseudoprocesses on the line with stable subordinators H^{t), t > 0, which possesses 
a Cauchy distributed law. As t oo the distribution (4.1) converges to the uniform 
law. 



Remark 4.3. In view of (4.2) we note that for 82„ yH^n (t) j ^ t > 0, the probability 
of staying in the right-hand side of has the remarkably simple form 



Pr{-|<e..«(,^A«,)<|} 



arctan e 



Vi > 0. 



(4.6) 
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Figure 3. In the picture the density of the circular Brownian 
motion (dotted Hne) and the kernel (4.1) are represented. 




We now pass to the Poisson kernel associated to odd-order pseudoprocesses. The 
asymmetry implies that the density of the composition Q271+1 

(^H^it)^ t > 0, 

is bit more complicated than (4.1). 



(4.7) 



-2a„t 



Theorem 4.4. The composition yH^"{i)j, t>0, has density 

Pr|e2„+i [H^it)) edej = —- 

and distribution function 

Pr{e2„+i (H5^(t)) <0} = 



2e~ 



' cos {0 + bnt) ' 



arctan tan - arctan I'^'^-Tt tan ^ 



arctan 



^ tan 



9+b„t 



n ^ e+b„t ^ ^ 



arctan 



[T-Zt tan^f, TT < 6* < 27r ^ . 

(4.8) 



whe 



cos 



bn = sin 



(4.9) 



2(2n+l)' "'"2(2n+l)' 
Proof. Let h_i_ (s, t), s, i > 0, be the density of a positively skewed stable process 
of order 2n+i ■ Then, in view of (2.2), we have that 



Pr |e2„+i (h^^ (t)) edej ^de + - ^ < 



{kO + k 



2n+l 



s) h_j_{s,t) 



= dd 

_ de 



1 1 

27r TT 



1 



^6^'^"*'= cos(fc(6' + 6„i)) 



k=l 
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2^r 



-2a„t 



2e- 



• cos [0 + bnt) ' 



In view of (4.4) we can write 
Pr{e2„+i <0} 

1 . 1 - 



2tt 



dy 



1 + e-2a,>t - 2e-""* cos (y + fe„t) 



arctan 



l+e' 
1-e- 



:;tan^ 



arctan 



-t tan 



6„t 



arctan 



l+e' 



"F tan 



9+fc„f 



- arctan nj-t^ 



tan 



2 ' 



(4.10) 



< 



27r 



(4.11) 

□ 




Figure 4. The distribution of the hitting point of a planar Brow- 
nian motion is obtained as subordinated circular pseudoprocess in 
the even case (Fig. 4a) and odd case (Fig. 4b). 



Remark 4.5. From (4.11) we arrive at the following fine expression 

Pr{e2„+i (i/5^(<)) <e} 

(l-e-2«^*)tanf (l + tan2%*) 

^2 taTi2 _ 



arctan 



(1 - e-«"*)^ + 4 tan f tan Y + (1 + e""-*)^ tan^ ^ 



0e [0,2^), 
(4.12) 

/rom which we are able to explicitely write for 82ji+i [H^"+^ (t) ) the probability of 
staying in the interval (0, tt) as 

Pr (o < e2„+i (t)) < vrl = - arctan ^ 

L \ / J vr sii 

w/ii/e for (O, ^) iwe obtain 

Pr{0<e2„+i (H5^(<)) < |} 



sinh Ont 
sin 6„ i ' 



yt > 0, (4.13) 
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2n+1 = 3 




2n+1 = 5 


\ / '• 

\ / 




\ / 








o 




6 Tt 2ti 
2n+1 =11 




n 2n 
2n+1 = 31 














o 





Figure 5. Distributions related to odd-order Poisson kernels (for t — \) 



arctan 



= — arctan 

TT 



(l+tan2 



(1 - e-<^r.tf + 4tan Y + (1 + e"''"*)^ tan^ 



sinh a„t 



2 sinh^ cos2 V + e'^" * sin 6„ < + 2 cosh^ sin^ V 



2 2 

1 sinha„i 

— arctan — — - — . 

TT cosh a„t — cos o„t + e""*^ sino„f 



(4.14) 



By means of the same manipulations leading to (4.14) we arrive at the alternative 
form of the distribution function for 9 G [0, tt], 



Pr |o < e2„+i (^H^ (t)) < ^} = ^ arctan 



sinha„t tan ! 



cosh ttnt — cos bnt + g""* sin bnt tan 1 



(4.15) 



Remark 4.6. In the third-order case we can arrive at the Poisson kernel (4.7) for 
n — \ by considering that (see [24]) 



Pr {x,[HHt))edx} _ r ds > X 



dx 



/3s 



3s/ sV3s 



/3s 



^/it 



27r 



L\ -X- 3*1 
2/ ^4 

(4.16) 



The wrapped up counterpart of (4.16) becomes for 9 G [0, 27r) 

oo 

Pr |e3 (i? 3 (t)) e ^6*} = Pr {Xg (t)^ e d{9 + 2m7r)} 

7n— — oo 

,n ^ r ds ..(9 + 2mn\ t ,.( t \ 
mtToo-'o v3s V v3s / sV3s Vv3s/ 

,y3i ^ t d9Vit °° '•^ 



:d9- 



2tt 



E 



m— — oo 



2m^+ |)' + 3f 



E 



m— — oo 



ds- 



(8 + 2,„„+|) _(V3t^ 

e 2s e 2s 



/27rs 



%/2^ 
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227r 



ds 



1 + 2^, 



fc=i 



k 



d9_ 
2^ 



1+2^:, 



cos 









k (^9 













e 2 



V27rs3 



27r 1 + e-v^* - 2e-^* cos (6* 



(4.17) 



5. Circular Fresnel pseudoprocess 
The equation of vibrations of rods 

^f{x,t) ^ ~K^^f{x,t), xeR,t>0, (5.1) 

is a sort of biquadratic heat equation which has suggested the construction of a 
pseudoprocess caUed Fresnel pseudoprocess (see Orsingher and D'Ovidio [23]) and 
denoted by F{t), t > 0. For the standard equation of vibrations of rods 

\f{x,0) = Six), 



the fundamental solution reads 



f(x,t) = ^=cos — - - . (5.3) 

For cylinders C — {x : aj < x{tj) < bj}, where x : t ^ x{t), t > 0, are the sample 
paths of F, we construct a sign-varying non-markovian measure /i based on (5.3) 
as (see formula 4.1 of [23]) 



m(c) 



''"^2.£^cos(vfc^ 



''3- 




The extension of the measure fi on the field generated by cylinders is performed in 
the same way as in the case of pseudoprocesses related to (1.1). We note that the 
composition F {\B{t)\), t > 0, has density coinciding with the fundamental solution 
to the fourth-order heat equation (Benachour et al. [6]). Furthermore F{Tt), t > 0, 
Tt — inf {s > : B{s) = t} has distribution coinciding with a bimodal Cauchy-type 
r.v. (see Orsingher and D'Ovidio [23]). 

We now study the counterpart on the ring, denoted by S{t), i > 0, of the Fresnel 
pseudoprocess F, and present its Fourier representation in the next theorem. 

Theorem 5.1. The density measure f{9, t) of the Fresnel pseudoprocess ^{t), t > 0, 
on the ring ^ is the wrapped up version of (5.3) 

,^ , 1 /(0 + 2m7r)2 7r\ , , 

^ m— — C30 

and has the Fourier expansion 

1 1 ^ 1.2+ 

f(0, t) = + _ y cos — cos ke. (5.6) 
27r TT ^ — ' 2 
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The densities (5.5) and (5.6) are solutions to 



at- 



^Ke,t) ^ -^^f{e,t), ee[o,2TT),t>o, 



[fie,o) = sie). 

Proof. The Fourier expansion of f{0,t) is 



f{9,t) = flo + ^flfe cos k6 



(5.7) 



(5.8) 



where 



2ir 



COS kO 



_j oo 

EGOS 



/2'Kt 



V27rf 

1 



E 



m+l 



COS 2fc7r0 COS 



{0 + 2m7r)2 TT 
2t 

/(27r6i)2 _ TT 



e**^" cos 



TT 

2t ^ 4 



ad = — cos . 

TT 2 



(5.9) 



This immediately leads to (5.6). The check that (5.5) and (5.6) are solutions to the 
equation of vibrations of rods (5.7) is straightforward. □ 

Remark 5.2. The law of the composition of the Fresnel pseudoprocess S(t), t > 0, 
with a positively-skewed stable process H^lt), t > 0, l3 G (0,1], has the following 
representation 



oo 

It) =^ + ^y 

' ' 27r 27r ^ 



cos k0 



fe=i 



1 1 
27r TT 



cos kOe 



y ^ cos — TT- sm — 



k=l 



V 2^ 2 



(5.10) 



Instead the composition of the Fresnel pseudoprocess d(t), t > 0, with the inverse 
L'^{t), t > 0, of the stable process H'^{t), t > 0, has the Fourier representation 



1 1 °° 

ZTT ZTT ^ — ' 



COS k6 



oo 



2~ 



- + -^cosfc^^;,., [-^j 



(5.11) 



(5.12) 



In (5.10) we used the characteristic function of positively skewed stable processes 
and in (5.11) we used (3.2). Each term of (5.11) is the solution of a time fractional 
diffusion equation while f'^{9,t) solves the fourth- order equation 



g2^ 

dt^ 



r{9,t) 



1 

'229^ 



no,t), 



0e [0,27r), t>0. 



(5.13) 



For v — ^ we extract from (5.12) the Fourier expansion of the fourth-order heat 
equation (1.13). For v = 1, since E2^i{—x) = cos^/x, we extract from (5.12) the 
fundamental solution (5.6) of (5.7). 
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5.1. Poisson kernels and circular EVesnel pseudoprocess. If we consider the 
Poisson kernel V (r, ^E*, i?, 0) where is a r.v. with Pr {5' — = ^ we have that 



1 
1 



2rR cos xjJ) 



27r (i?2 + - 2ri? cos cos ( 



(5.14) 



(2ri? sin-f/; sine 

The distribution (5.14) with respect to (/) can be interpreted as the distribution of 
the hitting point on the circumference of radius i? of a Brownian motion starting 
with probabihty ^ either from (r, +'0) or from (r, —?/')• The distribution (5.14) for 



^ = ±^ takes the fine form 



/ vr 



1 i?4 _ 

2^ i?4 + 7-4 + 2r2i?2 cos 20 ' 
For arbitrary values of ^' (5.14) can be written as 

1 — x—^ / 2r_R sin sin ( 



(5.15) 



¥.V {r,^,R,> 



2n i?2 _ ^2 _ 2r-i?cosi/' 



E 

fc=0 



i?2 -|- r2 — 2ri? cos "0 cos < 



2fc 



Remark 5.3. The composition 

din) , where n = inf {s > : B{s) = t} =™ H^V2t), t > 
can be viewed as 

din) ^^<B{zn)+^-^i-zn) 



(5.16) 
(5.17) 



2 

/2it 



1 



:«8 iJ 



/2it 



/2ii 



re. 



'2it 



(5.18) 



where 58 and O2 are circular Brownian motions differing only by a scaling factor. 
Thus by considering Theorem 4-1 we obtain for ^ (rt) 

j-V2t(l-i) ^ _ g-\/2t(l+i) 

^^^+'''> cos6i + e-\^(i+*)* 
(5.19) 



r{6,t) = 



1 

22^ 



1 

l-2e"75(i~^)cos6' + e-^(i-*)* ' 1 - 2e"75( 



Now we note that (5.19) presents the structure of the sum of two Poisson kernels 
written as 



p{9,t) = 



1 

2V 



(1 - rl) 



l-rl 



2ri cos t 



1 + To — 2r2 cos ( 



whe 



Ti — e ^ 



(5.20) 
(5.21) 
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After some manipulations we obtain 



p{0,t) 



2% 



8(1 - rir2)rir2 [1 + rir2 - (ri + coad] 



[22rir2 cos 6* - (ri + r2)(l + rir2)] - (ri - r2)^{l - rjrj) 



and by considering (5.21) we arrive at 
1 



(5.22) 



1 - e^^M n +e-^* - 2e ^^ cos 6*008^ 



2tt 



2e ^ cos - cos -7= ( 1 + e ^^ 



V2 



(^1 + e + (^1 - e-2\/2t^ sii 



.2 t 



72 

(5.23) 



We note that (5.22) for ri — — r2 = r coincides with (5.15) when R — I and 
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